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By starting with the Maxwell theory of electromagnetism, we study the change of polarization 
state of light transmitting through optically anisotropic media. The basic idea is to reduce the 
Maxwell equation to the Schrodinger like equation for two levels (or states) representing polarization. 
By using the quantum mechanical technique, the density matrix, and path integral, the evolution 
of the Stokes parameters results in the equation of motion for a pseudospin representing a point on 
the Poincare sphere. Two typical examples relevant to actual experiments are considered; the one 
gives the generalized Faraday effect, and the other realizes an optical analog of magnetic resonance. 

PACS number: 42.25.Ja, 78.20.Ek 



, The study of propagation of light (or an electromagnetic wave) in optical media has long been one of the major 
^ ' subjects in physics. We mention, for example, the classic monographs and the works on nonlinear optics ||Q. 
The characteristic quantity describing the light propagation is the concept of polarization. The study of polarization 
has also a long history |l],|2| , which forms a basis of modern crystal optics. The simple way to describe the polarization 
state is given by the Stokes parameters (or vector). Geometrically, the Stokes vector is realized as a point on the 
• so-called Poincare sphere. The Stokes vector or Poincare sphere play a powerful role for analyzing the change of the 
polarization state of light transmitting through anisotropic optical media As for the equation for evolution of 
the Stokes parameters, the phenomenological description has been known in the area of optics, which uses special 
C*~) | mathematical device such as the Jones vector or Miiller matrix 1^-0] . 

Having given a brief overview of the developments achieved so far, we address a novel formalism of the evolution 
for the polarization state of light transmitting through anisotropic media. Apart from the previous phenomenological 
approaches |^]-f?J , our theory is based on the first principle starting from the Maxwell theory of electromagnetism [Q , 
, where we use the more refined form than the original Maxwell equation. Namely, the Maxwell equation is reduced to 
the wave equation a la Schrodinger equation for two levels || , which is of first order in time (we call this the Maxwell- 
Schrodinger equation hereafter) and the dielectric tensor plays a role of Hamiltonian. By applying the technique used 
in usual quantum mechanics, such as density matrix as well as path integral, to the Maxwell-Schrodinger equation, we 
obtain the evolution equation for the Stokes parameters as an equation for a pseudospin which represents a point on 
the Poincare sphere. This is our main consequence. As typical applications of the equation of motion for pseudospin, 
we consider the polarization change in specific media, for which the dielectric tensor has the same structure as the 
O ■ Hamiltonian for a real spin in external magnetic field. Specifically we are concerned with two cases: The first example 
, is the pseudospin in uniform "magnetic field" which leads to the generalized Faraday effect. The second example 
is the pseudospin in oscillating as well as uniform field, by which we conjecture a possible occurrence of an optical 
rN | analog of the nuclear magnetic resonance (NMR) . These examples may be accessible to actual experiments and would 
enable us to reveal new aspects of polarization phenomena that have not been expected by previous works. 

Maxwell-Schrodinger equation. — We consider the plane electromagnetic (EM) wave of the wave vector k (k means 
the wave vector in the vacuum) travelling through the dielectric medium in the z direction. The medium is anisotropic 
with respect to the propagation direction and let e be the dielectric tensor. We assume z axis to be one of the principal 
axes of the dielectric tensor, namely, the direction along which the one of the eigenvalue of e. When the medium is 
isotropic, the eigenvalue is prescribed to take the value no. Thus the EM wave has the form like D(z, t) = D(z) exp[iwt], 
and the wave equation for the displacement vector D(z) is given by |jj 

where k = — . In the geometry under consideration, the dielectric tensor is taken to be 2 x 2 matrix. Under the most 
general condition that is governed by the external static electric and magnetic fields or mechanical constraint, e can 
be chosen to be a Hermitian matrix |Q, which means that the medium is transparent for the light transmission (we 
will consider elsewhere the case that there is an effect of absorption of light). Furthermore we consider the general 
situation that the medium is inhomogeneous, namely, the e depends on z. We now set the ansatz for the wave D 

D(z) = f (z) exp[ikn,Qz], (2) 
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where the amplitude f(z) is given by the 2x1 row vector, 



f = i (/i,/ 2 )=/ie 1 + / 2 e 2 . (3) 

f is a slowly varying function of z compared with the wave length, namely, k » |/'|, which implies that we consider 
the short range approximation. Here (ei, denotes the basis of linear polarization. In the short wave approximation, 
the amplitude is shown to satisfy the equation 

i\f z +{e-nl)f = 0, (4) 

where A = ^ is just the wave length in the medium of refractive index no divided by 2tt. Note that the second order 
differential term /" is discarded, since this is much smaller than the first order differential /'. In this way, the above 
equation can be regarded as an analog of the Schrodinger equation for two-level state, where A just plays a role of the 
Planck constant and z plays a role of the time variable. The components (/i,/2) of the vector f couple each other 
to give rise to the change of polarization and the "Hamiltonian" is given by h = e — n\. This form of h represents 
the deviation from the isotropic value, that is, "degree of anisotropy" , namely, the deviation governs the change of 
polarization state. From the hermiticity, the most general form of this is written as 



a ft + v) 
(3 — ij —a 



(5) 



Now for the later use, it is convenient to transform the basis of linear polarization f into the circular basis ||, that 
is, e± = -j= (ei ± ie 2 ), hence the Schrodinger equation becomes 

ixf z = 6*, (6) 



where ip = Ti = t (i/)*, ipZ), H = ThT 1 . Here T is the unitary transformation of the 2x2 matrix given by 

_L ( 1 * 

Thus the transformed Hamiltonian turns out to be 

H = ThT- 1 =( 7 . " + l/? V (8) 
V a - i(3 -7 J 

which is written in terms of the Pauli spin; H = y\ —1 hi<Ji. The formal solution of the above Schrodinger equation 
is given by ip(z) — T(z)ip(0) with T(z) being the evolution operator, 

f{z) = Pexp[-^ J H{z)dz\. (9) 

Here P denotes the path ordered product which is necessary to handle the z dependence of H . 

Density matrix and equation of motion of the Stokes parameters. — We now consider the reduction of the above 
Schrodinger equation. This is carried out by using the density matrix, for which we have two cases, the pure 
polarization and the mixed polarization (or partially polarized) state. Here we restrict ourselves to the former case 
in order to simplify the argument ptl| , thus the density matrix is defined as 

>=^=(M)' (10 » 

In terms of the components of the function t/j given above, we have the definition for the Stokes parameter [[TT|| ; 
Si = ?/^cri?/>, So = ip^lip, where u^i — 1,2,3) means the Pauli spin matrix. These variables satisfy the relation 
Sq = Sf + 5*1 + Sf, which is equivalent to the equation detp = 0. Furthermore, in the case that the Hamiltonian is 
Hermitian, we can adopt the conservation of probability ip*ipi +^2^2 = 1- S° if we use the spinor parametrization 

9 9 

?Ai=cos-, ip2 = sia.-exp[i<f>], (11) 
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we have p = i(l + So - ), where the vector S = (S x , S y , S z ) is given by 

S x = sin cos 0, S y = sin sin 0, S' 2 =cos#. (12) 

This forms the Stokes vector and is described by the point on the Poincare sphere. We illustrate some typical values: 
(i) 9 — 0; the north pole that corresponds to the left handed circular polarization, (ii) 9 = ir; the south pole that 
corresponds to the right handed circular polarization, (iii) 9 = ^; the equator which represents the linear polarization. 
The equation of motion for the density matrix is written as 

i\f z = [H,p}. (13) 

Here [H,p] = Hp — pH. Using the commutation relation for the Pauli spin, [<7i,<7j] = 2ieijkO~k, we can deduce the 
equation of motion for the pseudospin from ( ft3| ) 

A? = SxG, (14) 
dz 

where the effective "magnetic field" is defined as G = (2a, 2/3, 27). If we introduce the "classical" counterpart of the 
Hamiltonian (||) 

H = 2aS x + 2(3S y + 2-yS z 

= 2a sin 6 cos (j> + 2(3 sin 9 sin <j> + 2^ cos 9, (15) 

we have an alternative form of the equation of motion p2| 

* 1 OH ; I OH 



XsiaO d(j) Asinf? 89 

The equation of motion (|I7J) can also be obtained as a result of the asymptotic limit of A — * |l3| . This may be 
achieved by the fact that the set of states for pseudospin, ip = \ip > forms a Bloch state that satisfies the completeness 
relation: J \ip > dfJ, < tp\ = 1 with measure dp, — sm9d9d(f> (just the volume on the Poincare sphere). Let us consider 
the transition amplitude that is given by sandwiching the evolution operator (^|) with two initial and final spin states. 
By adopting the procedure of "time slicing" and inserting the completeness relation at each time division, we get the 
path integral expression 113] 

< V/|f(z)|lfc >= J exp jis[V*,V]}£>M(V>) (17) 
with the path measure Dp(ip) = Y[ z dp[ip*(z), iji(z)] and S is the "action function" 

S = I < ip\i\4- ~ H\tb >dz. (18) 
J dz 

In the limit of A — > 0, we have the stationary phase condition SS = leading to the equation of motion for the 
pseudospin, i.e., (16). 

Typical applications. — We shall consider some special cases that can be described by the general formalism, (i) 
We first consider the model for which the dielectric tensor depends on the external magnetic field as well as electric 
field. The kinematical symmetry implies that h has the form 



a ij 
-17 —a 



(19) 



Here 7 is proportional to the uniform magnetic field (the strength is B) applied in the z direction; 7 = gB ]lq ]. Thus 
according to the formula (||), it is transformed to H = ^a z + aa x . This is further transformed to the Hamiltonian for 
the spin in uniform field of strength T, that is applied in the z' direction by rotating about the y axis by an amount 
of the angle i], such that 

S x = cosi]S x — smr]S z , S' z = sinr)S x + cos i]S z , (20) 
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where V — \/"i 2 + a 2 together with the angle cos 77 = 7/r. Thus the equation of motion for the pseudospin becomes 



dS' x _2T dS y _ 2T 

dz ~ A by ' dz ~ A 1 j 

for which we get the solution S' x — sin(^-z + 6q), S y — cos(^z + Oq) and S' z = 0. In terms of the original spin, it 
gives 

2r 2r 

S x = cosr?sin(— z + 6q), S z = — sin 77 sin(— z + 9 ). (22) 
A A 

This means the following feature, if the light wave is initially in the linear polarization, then it turns out to be an 
elliptic polarization after transmitting through the medium (see Fig. Q). As a special case of a = 0, the initial 
polarization changes to S x — sin(^z + 9o), S y — cos(^z + 9q), namely, the initial polarization plane simply changes 
by an amount of the angle (j) = -?-l after travelling the distance z — I. This rotation of the polarization plane is known 
to be the Faraday effect [see Fig. ^|(a)]. 

(ii) We consider another example for which the polarized light propagates in the medium such that the dielectric 
tensor has a periodical structure besides the effect of the external magnetic field as in the case (i). For such a system 
the tensor h is given by 

h=( 7 ° coswz 7 ° sinwz )+( (23) 

V7oSino;z — -focosuz J V — 17 I 

which is transformed to H(z) = G(z)»S, where the pseudomagnetic field has the component: 

G(z) = (270COSWZ, 27osinwz, 27), (24) 

namely, a static field along the z axis plus an oscillating field rotating perpendicular to it with the frequency to. This 
feature is familiar in magnetic resonance. Using the classical counterpart of the above Hamiltonian, which is given as 
H(z) = G(z) ■ S = 270 sin 9 cos(</> — ujz) + 27 cos 9, the equation of motion is derived by using the general formula jl^ ) 

a • (a \ 

a = — sm(0 — ujz), 

A 

2 

<j> = — -r [70 cot 9cos(<l> — loz) — 7]. (25) 
A 

One sees that this form of equations of motion allows a special solution 

(j) = ujz, 9 = 9q(= const), (26) 
where the following relation should hold among the parameters #o 5 7i7o : 

cot* =(^-^). (27) 
\7o 270/ 

Equation (|26| ) may be called the "resonance" solution, since it corresponds to the solution for the forced oscillator. 
The set of parameters (7,70,0;) satisfying (E7|) for a fixed value 9 — 6q belong to a family of resonance solutions. 
Indeed, this set of parameters forms a surface in the parameter space (7,70,0;), which we call the "invariant surface" 
and characterizes the resonance condition. The condition ( p7j ) is crucial, since all the quantities on the right hand 
side are given in terms of constants that may be allowed to be compared with experiment. The physical meaning of 
the above invariant surface is as follows: If given is an initial elliptic polarized wave with the angle 9q satisfying the 
resonance condition, then there is no change in its shape during its transmission. Only its axis rotates by an amount 
of loI that is governed by the period inherent in the dielectric tensor [see Fig. 0(b)]. 

Here we give a remark on possible realization of the periodic structure of the dielectric tensor in actual systems. 



The one realization is made by using special materials, such as cholesteric liquid crystals |16|. The periodic structure 
is naturally realized by the helix inherent in liquid crystal. Indeed, the dielectric tensor can be given by the form of 
the first term of (|23) p!q |. Another realization may be given by using the mechanical one, namely, let us consider 
the elastic body under the pressure that is periodically modulated, then this causes the periodic modulation of the 
dielectric tensor according to the procedure known in the mechanical-optical effect [[jj. The periodical oscillation of 
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the pressure may be generated, for example, by the piezoelectric effect. Now having settled the periodic structure, 
we take into account the term coming from the uniform magnetic field that is applied along the same direction as 
the axis of helix. We are reminded of the analogy with the NMR; if the initial state is in the left handed circular 
polarization, which corresponds to the state of spin-down, the probability for transition to the state of spin-up (right 
handed circular polarization) is given by 

P d - r) = ( T 2 ° } !, M2 sin 2 \ Az, (28) 

(si - uj) a + (j^y 2 

where 7q = 3£ and SI = % and A 2 = (SI — u>) 2 + (jq) 2 - If the magnetic field is chosen so as to synchronize the period 
of the helix, we expect an analogous effect with the magnetic resonance. 

This work was inspired by the discussion at seminar class that had been instructed by one of the authors (H.K.). 
The authors thank the attendees of the seminar. 
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FIG. 1. The initial linear polarization, which is marked FIG. 2. (a) The change of linear polarization by an 

by /, changes to the final elliptic polarization are marked amount of cj) = %t-l, (b) The change of the axis of elliptic 

by F. polarization by an amount of u)l. 
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